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Abstract The relationship between the change of forest resources and climatic factor in the “Three-North” re¬ 
gion of China were studied in this paper. The predicting equations of climatic factor (dependent variable) with re¬ 
gional independent variable (longitude, latitude and altitude) and stand independent variable (forest coverage rate), 
were developed by extensively using the linear and nonlinear regression methods. With these models, we can cal¬ 
culate the ecological benefit of Shelter-belt forest. 
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Introduction 

In order to solve the serious ecological problem we 
faced and improve the phenomenon of land desertifi¬ 
cation in the north area of China. China has begun to 
construct the Three-North Protective Forest System 
since 1978. With the hard work of twenty years, it is 
told “The greatest ecological project of the world” and 
“The green great wall”. Although the country invested 
lots of capital to build the Three-North shelter-forest 
system through many ways such as Three-North 
special investment, finance appropriate funds, finan¬ 
cial loans and take out funds for bring up trees, but 
lack of funds was a very seriously problem during the 
process. So it is necessary to levy expenses of bring 
up ecological resources from the units and depart¬ 
ments which benefit from the protective forest, and 
give some compensate to the ecological benefits 
bring from the protective forest. In this paper, the 
relationship between the changes of forest resource 
and climatic factors in Three-North district was stud¬ 
ied, and the regression equations of climatic factors 
and regional factors were developed. The results can 
provide some scientific bases to calculate the 
ecological benefits of protective forest. 

Materials and methods 

Data collection 

The data used in this paper was collected from more 


than one hundred counties of Heilongjiang, Jilin, 
Liaoning and Inner Mongolia. These areas lay be¬ 
tween longitude 122°12’~130°18’ E and latitude 
43°~50°24’ N. The data contains geographic data 
(longitude, latitude, altitude) and forest coverage rate 
and climatic factors (mean annual temperature, an¬ 
nual precipitation, mean annual wind speed, annual 
evaporation, annual accumulated temperature for 2 s 
10 °C) of typical years of each county. 

Variables 

Three types of variables were used. 

Dependent variables: 

Y,= mean annual temperature, °C 

Y z = annual precipitation, mm 

Y 3 = mean annual wind speed, m/s 

Y„= annual evaporation, mm 

Y 5 = annual accumulated temperatureSMCTC, °C 

Regional independent variables 

X,= longitude ( *) 

X 2 = latitude ( *) 

X 3 = altitude, m 

Stand independent variable 

X 4 = forest coverage rate, % 

The mean, maximum, and minimum values of each 
variable were listed in Table 1 

Models 

In order to develop the models to express the regres¬ 
sion relationship of these variables perfectly, the 
following methods were used to build the regression 
models of each variable. 
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Table 1. The characteristic values of each variables 

Variables 

Mean 

Minimum 

Maximum 

Longitude ( ’) 

125.3620 

122.200 

130.300 

Latitude ( °) 

46.0970 

43.000 

50.400 

Forest coverage rate, % 

16.5640 

3.400 

71.100 

Altitude, m 

220.5650 

81.200 

528.500 

Mean annual temperature, 'C 

3.4663 

-0.500 

6.000 

Annual precipitation, mm 

495.6278 

344.200 

889.000 

Mean annual wind speed, m/s 

3.7543 

1.800 

5.200 

Annual evaporation, mm 

1507.5420 

752.500 

2083.300 

Annual accumulated tern- 

2723.7480 

2153.000 

3138.800 


perature S=10’C, *C _ 

We choose seven models and calculate the relati¬ 
on coefficient of every dependent variable and each 
independent variable, and according to these coeffi¬ 
cients to make sure of the relationship of dependent 
variables and independent variables. The calculation 
results were listed in Table 2 


According to the relationship of variables we can 
build regression models of each dependent variable 
and every independent variable as follows: 

Mean annual temperature: 

Y,=a+bX,+cX 1 2 +dX 2 +eX 3 +fX 3 2 +gX 4 + hX 4 2 +iX 4 3 
Annual precipitation: 

Y 2 =a+bX, +cX, 2 +dX 2 +eX 2 2 +fX 2 3 +gX 3 +hX 3 2 +iX 4 +jX 4 2 
Mean annual wind speed: 

Y 3 =a+bX 1 +cX 1 2 +dX 1 3 +eX 2 +fX 2 2 +gX 2 3 +hX 3 +l*exp(jX 4 ) 
Annual evaporation: 

Y 4 =a+bX.,+cX 1 2 +dX 1 3 +eX 2 +fX 2 2 +gX 3 + hX 3 2 +iX 4 +jX 4 
Annual accumulated temperature for ^10”C: 
Y 5 =a+bX 1 +cX 1 2 +dX 1 3 +e*exp(fX 2 )+gX 3 +hX 3 2 +iX 3 3 + 
J*exp(kX 4 ) 

Where: a, b, c, d, e, f, g, h, i, j and k are the parame¬ 
ters to be estimated. 

Estimate and test the parameters in these models, 
the result is list in Table 3. 


Table 2. Regression equation and relation coefficient of every dependent variables and each independent variables 


Dependent variables | 


Mean annual temperature 


Independent variables ( 

Longitude 

Latitude 

Altitude 

Forest coverage rate 
Longitude 
Latitude 
Altitude 

Forest coverage rate 

Longitude 

Latitude 

Altitude 

Forest coverage rate 
Longitude 
Latitude 
Altitude 

Forest coverage rate 

Longitude 

Latitude 

Altitude 

Forest coverage r ate 


Regression equations 
Y=a+bx+cx 2 
Y=a+bx 
Y=a+bx+cx 2 
Y=a+bx+cx 2 +dx 3 
Y=a+hx+cx 2 
Y=a+bx+cx 2 +dx 3 
Y=a+bx+cx 2 
Y=a+bx+cx 2 
Y=a+bx+cx 2 +dx 3 
Y=a+bx+cx 2 +dx 3 
Y=a+bx 
Y-a+b*exp(cx) 
Y=a+bx+cx 2 +dx 3 
Y=a+bx+cx 2 
Y=a+bx+cx 2 
Y=a+bx+cx 2 
Y=a+bx+cx 2 +dx 3 
Y=a+b*exp(cx) 
Y=a+bx+cx 2 +dx 3 
Y=a+b*exp(cx) 


Relation coefficient 
0.608276 
0.956242 
0.341027 
0.7 00556 
0.641864 
0.591687 
0.652937 
0.685769 
0.741650 
0.727203 
0.441814 
0.770203 
0.530595 
0.588132 
0.473286 
0.637882 
0.649570 
0.906724 
0.518169 
0.725672 


Annual precipitation 


Mean annual wind speed 


Annual evaporation 


Annual accumulated temperatureSMO'C 


Table 3. Estimated value of parameters 


parameters _ 


Dependent variables 
Y, 








Journal of Forestry Research, Vol. 10, No. 1, Mar. 1999 


43 


Form these statistical results, we can find that the 
regression relationship of these independent vari¬ 
ables were very significantly. So we can develop 
regression models as follows: 

(1) Mean annual temperature (V,) : 

V, =965.3284-14.4886X,+0.056858X, 2 - 

0.810199X 2 + 0.010425X 3 -0.000201X 3 2 - 
0.020214X 4 +0.00006X 4 2 

(2) Annual precipitation (V 2 ) : 
y 2 =-10207.1+808.2876X 1 -3.12181X 1 2 -2655.69X 2 + 

56.50568X 2 2 - 0.402962X/+1,642257X 3 - 

0.002303X 3 2 +2.863701*374.789X 4 - 

0.035938X 4 2 

(3) Mean annual wind speed (Y 3 ) ■. 
V 3 =-34632.7+889.8219X,-7.01699X, 2 +0.018437X, 3 - 

192.044X 2 +4.148953X 2 2 -0.029846X 2 3 - 
0.016956X 3 + 42.1575*exp(-2.19089X 4 ) 

(4) Annual evaporation (y 4 ) : 
y 4 =3055230-72197.5X 1 +568.5245X 1 2 -1.49322X 1 3 + 

252.9407X 2 -3.37242X 2 2 -0.943575X 3 + 

0.001048X 3 2 -11.5772X 4 +0.140510X 4 2 

(5) Annual accumulated temperatureS=10”C (Y s ) : 
Y s = 1028.517-608.601X,+10.55819X, 2 - 

0.04413X, 3 -0.001737*exp(0.26295X a )- 
6.7181 6X3+O.O21 269X 3 2 -0.00002X 3 3 - 
0.000352*exp (0.141712*X 4 ) 

The relationship between observed value and es¬ 
timated value of each dependent variables can be 
expressed by Fig.1 ~5. 



Predicted Values 


Fig. 1. Observed value and predicted value relationship 
of mean annual temperature 

Conclusion 

According to these models we can find that: 

(1) The relationships between the mean annual 
temperature and longitude, latitude, altitude or forest 
coverage rate are negative correlation, but the affec¬ 
tion of longitude, latitude or altitude to mean annual 


temperature is larger than that of forest coverage 
rate. 



Fig. 2. Observed value and predicted value relationship 
of annual precipitation 



Fig. 3. Observed value and predicted value relationship 
of mean annual wind speed 



Predicted Values 

Fig. 4. Observed value and predicted value relationship 
of annual evaporation 
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Predicted Values 

Fig. 5. observed value and predicted value relationship 
of annual accumulated temperature >10'C 

(2) The relationships between annual precipitation 
and longitude, altitude or forest coverage rate are 
positive correlation. The relationship between the 
annual precipitation and latitude is negative correla¬ 
tion. These four factors have obviously affection to 
annual precipitation. 

(3) The mean annual wind speed and these four 
factors are all negative correlation. The affection of 
forest coverage rate to mean annual wind speed is 
the largest. 

(4) The annual evaporation and longitude or lati¬ 
tude are positive correlation. The annual evaporation 
has negative correlation with altitude or forest cover¬ 
age rate. The affection of forest coverage rate to 
annual amount of evaporation is the largest. 

(5) The relationships between annual ^10"C ac¬ 
cumulated temperature and these four factors are 
negative correlation. The affection of the first three 
factors is larger than that of forest coverage rate. 

From the relationships of factors we can conclude 
that: 

The increase of forest coverage rate can cut down 
the mean annual wind speed of planting area, reduce 
the damage of wind and sand, and benefit to agricul¬ 
ture. Meanwhile, since the increasing of forest cover¬ 
age rate can cause the increases of the annual pre¬ 
cipitation and annual evaporation, it will cause little 
affection to water balance. 

The increase of forest coverage rate may cause 


the decreases of mean annual temperature and an¬ 
nual accumulated temperature for^10”C, which may 
cause a little negative effect on agriculture. 

As what is mentioned above, the change of stand 
factor (forest coverage rate) affects the regional cli¬ 
matic factors very obviously. The regression equa¬ 
tions can be set up depending on the relationships of 
the regional independent variables, stand indepen¬ 
dent variables and the climatic factors (dependent 
variables). Using these models we can calculate the 
ecological benefit of shelter-belt forests. 
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